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Abstract 

We study q-stars with one or two scalar fields, non-abelian, and 
fermion-scalar q-stars in 2 + 1 dimensions in an anti de Sitter or flat 
spacetime. We fully investigate their properties, such as mass, particle 
number, radius, numerically, and focus on the matter of their stability 
against decay to free particles and gravitational collapse. We also 
provide analytical solutions in the case of flat spacetime and other 
special cases. 
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1 Introduction 



The investigation of scalar field configurations coupled to gravity started 
with the work by Kaup, jT], and Ruffini and Bonazzola, [2]. They regarded a 
complex scalar field with no self-interactions. Quartic self interactions were 
taken into account in other works EJ. Scalar fields with a local U(l) 
symmetry were also investigated jHj- The results were generalized in a series 
of papers in other gravity theories [UJ |Hl El • 

Soliton stars appeared in the literature as stable field configurations con- 
sisted of abelian or non-abelian scalar fields, containing sometimes fermions 
or gauge bosons coupled to gravity [TUl HH UHl [201 IZH 123 123] - These 
configurations are stable even in the absence of gravity. A certain class of 
soliton stars are very large field configurations with radius of the order of 
lightyears fUl El Ej • They consist of a very large interior with very low en- 
ergy density of "kinetic" type resulting from the time variation of the scalar 
field. On the other hand the potential energy within the soliton is considered 
to be negligible. Also, the metric and matter fields vary very slowly in the 
soliton interior. A considerable amount of energy, comparable to that stored 
in the huge interior, is contained in the thin surface. Within the surface both 
the potential and the energy resulting from the radial variation of the matter 
field contribute to the total surface energy, but the contribution of the time 
variation of the field is very small, as it is supposed to rotate in its internal 
£7(1) space extremely slowly. 

Another class of non-topological soliton stars, the q-stars, appeared as 
relativistic generalizations of a certain family of non-topological solitons, 
namely q-balls. Q-balls have a certain role in particle theory and especially 
in the baryogenesis through the flat MSSM directions [TZj. They appear in 
Lagrangians, [EH d EEH1 EH] , with a global 17(1) or SU{3), SO (3) symme- 
try, or a local U(l) [TS|, when the scalar field takes the special value which 
minimizes the U / '(f) 2 quantity. 

The properties of the q-stars have been studied thoroughly in ^JJ El 
122 123 E]- There are q-stars with one or two scalar fields, [THl 122] , q-stars 
with non-abelian fields, [22], and q-stars with a fermion and a scalar field, 
[213 El 122] • These objects contain a large interior within which both metric 
and matter fields vary smoothly. The energy within the interior is the sum 
of the potential energy and of the "kinetic" energy resulting from the time 
variation of the matter field. In the case of q-balls and, consequently, the 
case q-stars the potential is always positive. We will name U, W and V 
the potential energy, the energy resulting from the time variation and the 
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energy resulting from the space variation respectively. There is a crucial 
relation that roughly describes the q-solitons and differentiates them from 
other classes of non topological solitons, namely: 



where u is the frequency with which the soliton rotates in its internal U(l) 
space, or, with proper generalizations that we describe, in more complicated 
spaces, is a typical value of the or one of the scalar field(s) within the soliton 
and m is the mass of the free particles. This relation means that W is of 
the same order of magnitude with U, in contrast with other non-topological 
solitons for which the U is negligible in the interior and W negligible within 
the surface. The common feature between q-stars and other soliton stars is 
the smooth variation of the matter field in the interior and the rapid change, 
within the surface, from a <fto value to zero. In non-topological soliton stars, 
the energy density within the surface is huge compared to the energy density 
in the interior, and so, although the surface is very thin, the energy, contained 
within it, is comparable to the energy stored in the interior. This does 
not hold true for the q-stars. They have an approximately constant energy 
density in the interior, the sum of W and U. At the surface, these quantities 
retain the same order of magnitude and V in now added with the same order 
of magnitude as well. Because the energy density is everywhere of the same 
order of magnitude, either in the interior, or in the surface, and because the 
surface is of width of m -1 , the total energy of the thin surface is negligible. 
It has been proved that q-stars are smaller than non-topological soliton stars 
but larger than boson stars with no soliton features in the absence of gravity. 

Gravity in 2 + 1, [2JJJ ED], dimensions coupled to scalar fields is studied 
for many reasons. One of them is the theoretical interest for a theory qual- 
itatively different from the corresponding theory in 3 + 1 dimensions. The 
solutions in a 3— dimensional theory of gravity are much more simpler be- 
cause there is no A/r 2 — 1 term in the Einstein tensor, where A is a metric 
field as we will see. This can lead even to analytical solutions for special 
field configurations composed of scalar fields (or fermion-scalar) coupled to 
3— dimensional gravity [21>1 1271 128]. 

The aim of the present article is to study the formation of q-stars in 2 + 1 
dimensions in anti anti de Sitter or flat spacetime, as a limiting case. We 
want to investigate the matter of their stability with respect to gravitational 
collapse and to fission into free particles and to study the influence of the 
global spacetime curvature in the features of the q-star. We study four dif- 
ferent kinds of q-stars, namely, with one and two scalar fields, non-abelian 
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q-stars and fermion- scalar q-stars. We thoroughly investigate their proper- 
ties, such as soliton radius, mass, particle number, as functions of either the 
cosmo logical constant, which is an absolutely independent parameter, or as 
functions of their internal frequency, a quantity with crucial role in the theory 
of non-topological solitons. We also give analytical solutions for the impor- 
tant case that the cosmological constant tends to zero and also for other 
cases, which show a very close similarity in their behavior, when compared 
to solutions obtained with the usual numerical methods. 



2 A q-star with one scalar field 

We consider a static, spherically symmetric metric: 

ds 2 = -e v dt 2 + e x dp 2 + p 2 da 2 , 



(2.1) 



with g tt = —e u . We regard a scalar field with the minimum-energy time 
dependence: 

<j>(p,t) = a(p)e-^ . (2.2) 

The action in natural units for a scalar field coupled to gravity in 2+1 di- 
mensions is: 

~R - 2A 



167rG 



+ g^(d„<f>)*(d u( f>)-U 



(2.3) 



where A stands for the cosmological constant, regarded here to be negative, 
or zero, as a limiting case. The energy-momentum tensor is: 

Tfj,v = {d,<t>)\dA) + {dMdAY - g^[g aP {d a <t>)*{d p <t>)\ - g,M . (2.4) 

The theory need not be fundamental and, thus, renormalizable but effective. 
The Euler-Lagrange equation for the matter field is: 



dU 



w\9Mv\9\<rd v ) - 2 

dm 



taking now the form: 



= , 



dU 



a" + [1/p + (1/2) {y' - A')K + eWe^a - e x ^-a = 



da 2 



(2.5) 



(2.6) 



The Einstein equations are: 



(2.7) 



With the above assumptions the Einstein equations take the form: 

y 

— e~ A — = 8tvG(—W — V — U) — A , (2.8) 
Ip 

v' 

e~ x — = 8tiG(W + V — U) — A , (2.9) 



where 



2p 



\ ■ 1 — 1 ( -V.' 2 . 



dp J \dp 



(2.10) 



There is a Noether current: 

= V^gTii^d^ - 09,0*) . (2.11) 

The current is conserved according to the equation: 

3% = • (2.12) 

The total charge is defined as: 

Q = J d 2 xj° , (2.13) 

now taking the form 

Q = 4tt / pdpua 2 e- u/2 e x/2 . (2.14) 



This charge represents the total particle number if the charge of each particle 
is unity. In the present case the energy of the free particles with the same 
charge is mQ, where m is the mass of the free particles. We rescale all the 
parameters of the Lagrangian with respect to the mass, and set m eqaul 
to unity. So, if the total charge is less than the total mass of the field 
configuration, the star is stable with respect to decay into free particles. We 
define for convenience: 

B = e- U , A = e~ x , (2.15) 
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and rescale: 



p = pm, uj = uj/m, a = a/m 1/2 , 
U = U/m 3 , W = W/m 3 , V = V/m 3 , 

where m in generally is a mass scale, here the mass of the scalar field. We 
also make the rescalings 

A87rGm 3 = A, f = V8nGmp . (2.17) 

Inside the soliton both metric and matter fields vary very slowly with 
respect to the radius. This leads to a considerable simplification to both Ein- 
stein and Lagrange equations. We see that V oc (f)(0) 2 / P 2 oc m 3 cx(0) 2 / P 2 oc 
Gm^a 2 / R 2 which is very small compared to the other energy quantities. P 
is the (unrescaled) soliton radius and P and R correspond to the rescaled 
radii, according to the rescalings under |2~T^1 and I2~T7I respectively 

So, ignoring the 0(Gm 4 ) quantities and dropping the tildes the Einstein 
equations take the form: 

~~T~~~ = u 2 o 2 B + U + A , (2.18) 
dr 2r 

_AdBl_ = _ u _ A 

B dr 2r 

and the Euler-Lagrange equation is: 

2 da 

A rescaled potential admitting q-ball type solutions in the absence of gravity 
is: 

U = |0| 2 - |0| 4 + ^|0| 6 = o 2 (\ - a 2 + . (2.21) 

The solution to the eq. 12.201 is: 

a 2 = l + cuVB (2.22) 

and U is: 

U = -(1 + lu 3 B 3 / 2 ) . (2.23) 
3 
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Within the thin surface the metric and matter fields change rapidly. We 
can write the Lagrange equation in the following form (unrescaled quantities): 



A 1 A dB\ do 1 5(U — W — V) 
~p~ ~ 2B^lp~ J ~dp ~ 2 To ' ^' ' 

There is a certain form for the gn{r) metric component (unrescaled quanti- 
ties), [SI], |22!: 

-fi'W = Mp) = 1 - ^ - (g _ 2 2 ^_ 1) , (2-25) 

where _E P is the total energy density contained within a sphere of radius p and 
D is the spacetime dimensionality, here equal to 3. The above relation can 
be verified numerically and will be discussed more extensively when dealing 
with q-stars with two scalar fields, which offer a richer variety of analyti- 
cal solutions. Using the relation holding within the surface approximately 
namely the: A ~ 1/B and rescaling in the usual way, we find that: 

IS(U-W-V) j 1 - - £ d?re - ±f£ \1 d 2 re - 2Af 2 , 

— 2 '- ~ V87iGm 2^2 g drJo m 4 , 

2 bo r 

(2.26) 

where we regarded that within the surface the field decreases from a Oq value 
to zero and that oq ~ m 1 / 2 . Also, the surface width is ~ mT 1 . So, dropping 
the 0{y/Gm) quantities, we find: 

s { w-u-v) =0 

do 

Because at the outer edge of the surface all energy quantities are zero, we 
find a first integral to the above equation: 

V + U-W = 0. (2.28) 

Eqs. 12. 26112.28*1 hold true only within the surface. At the inner edge of the 
surface o' is zero in order to match the interior with the surface solution. So, 
at the inner edge of the surface the equality U = W, together with eqs. 12.221 



I2.23| gives the eigenvalue equation for uo: 

0-1/2 ,1/2 
= = ff«L . (2.29) 

2 2 y ! 

which has the right limiting value when gravity is absent, i.e.: A sur — > 1. 
So, uj, apart from a certain parameter referring to the soliton properties, is 
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also a measure for the gravity strength, equivalent to the metrics. The same 
discussion, eqs. I2.24H2.21JI holds true in any case of q-star. 

Solving the Einstein equations and in order to find the soliton parameters, 
we use the G° component of the Einstein tensor for the calculation of the 
energy or, equivalently, the form 

E = 7r(l — A(r) - Ar 2 ) , r -> oo , 

(rescaled quantities) arising from eq. 12.251 and eq. 12.141 for the calculation 
of the charge. The thin surface has a negligible energy and charge contri- 
bution. The independent variables are two: The cosmological constant, a 
global feature of the spacetime, and the soliton frequency. 

Outside the soliton the energy quantities are all zero. Einstein equations 
can be solved analytically: 

A(r) = A sm - Ar 2 + AR 2 = 1 Ar 2 , (2.30) 

7T 

, . A SUI - Ar 2 + AR 2 ( E , 9 \ 1 . . 

B( r ) = — — — = 11 Ar 2 , (2.31) 
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sur \ 

with R the star radius and E the total mass. The same relations hold true 
for every case of q-star we discuss. 

We will now give the analytic solution holding true when A — > 0. In this 
case the metric B is constant everywhere and the other parameters of the 
soliton are given in the following equations: 

B(r) = 1/A sm = JL , (2.32) 
A(r) = 1 - \r 2 , (2.33) 



R= ) J±(l-A mt ) = ^(l-4 U *) , (2.34) 
E = 7r(l — A sur ) = 7r(l — Auj 2 ) , (2.35) 

Q = 4vr(l - y/A^) = 4tt(1 - 2u) . (2.36) 

In the above equations R is the soliton radius determined when A(r) = 
l/B(r) and A sm is the value of the A(r) at the surface of the star. There 
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2 4 6 8 10 12 

Charge 

Figure 1: Energy as a function of charge for several values of the cosmo logical 
constant for a q-star with one scalar field. The numbers within the figures 
HJ IHl EU and IT7II2UI denote the absolute value of the cosmo logical constant. 
We start from a large value of the metric at the surface, A sur , near unity, 
corresponding to the weak gravity limit . The values of the energy and 
charge are small. As we decrease the metric the values of energy and charge 
increase. If A ^ and when reaching a certain value of A sur , then energy and 
charge decrease. The initial increase in the energy and charge values is not 
visible in the A = —2.5 case, because it takes place nearly to the maximum 
energy and for a short range of the surface gravity parameter. 
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1.2 



-10 



Cosm. constant 



Figure 2: The matter field value at the center of the soliton as a function of 
the cosmological constant for a q-star with one scalar field for four different 
values of -A sur . The frequency can be obtained using relation 12.291 For zero 
cosmological constant the field is everywhere constant within the soliton, 

1 /2 

equal to (1 + uj/Allr) = 1.5 1 / 2 . As the cosmological constant increases in 
absolute values the field takes larger values as a reflection of the increase 
in the gravity strength. Also, the field has larger values when u decreases, 
because small uj indicates stronger gravity. 




Figure 3: The soliton radius as a function of the cosmological constant for a 
q-star with one scalar field for four different values of A sur . 
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1.5 

: 



-10 -8 -6 -4 -2 

Cosm. constant 

Figure 4: The total soliton mass as a function of the cosmological constant 
for a q-star with one scalar field for four different values of A snr . 
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Figure 5: The particle number as a function of the cosmological constant for 
a q-star with one scalar field for four different values of A sur . We find that ra- 
dius, energy and charge increase for small absolute values of the cosmological 
constant but for larger values they decrease. 
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are some interesting results obtained when investigating the properties of the 
analytical solution. The total soliton energy is: 

E = —R 2 . (2.37) 
4 

This agrees with the intuitive result that the total energy is analogous to the 
"volume" of the field configuration. Also, E/Q < 1 = the rescaled mass of 
the free particles for every < A sur < 1 for the case of zero cosmological 
constant. This means that there is no possibility for the soliton fission to 
free particles, because this procedure is energetically forbidden. From the 
a 2 = 1 + ujB 1 / 2 relation and eq. 12.321 we find that in a flat spacetime, the 
scalar field has an everywhere constant absolute value, here equal to 1.5. 
So, when A sm — > 0, i.e.: when a horizon is going to be formed, there is no 
anomaly at the center of the star. 

In our figures, every dashed line corresponds to an analytical solution. 
The total energy is in 8ttG units, the particle number in 8irGm units, with 
rh — 1 and the soliton radius in (SirGrh 3 ) 1 ^ 2 units. In figures |2E starting 
with an initial zero value of the cosmological constant, we increase its absolute 
value and terminate our calculations when the soliton energy equates to the 
energy of free particles with the same charge. In figure Q we start from a 
large value of the metric field A SVLI , near unity, or, equivalently, a large value 
of the characteristic soliton frequency, near 1/2 (which according to eq. 12.291 
is the maximum limiting value of the frequency when gravity is absent, i.e., 
when A sur — > 1) and gradually decrease the frequency. 

We can observe, as a general comment of our numerical results depicted 
in figures [TJIH1 that the effect of the increase, in absolute values, of the cosmo- 
logical constant causes a consequent increase in the energy, charge and radius 
of the q-star. This is a counter-balancing effect: Negative cosmological con- 
stant means that different particles within the soliton, corresponding to the 
4> field, tend to move along deviating geodesies. So, the increase in mass 
and the other soliton parameters is necessary in order the object to avoid 
this kind of "deviation", or, in other words, to generate "positive" gravity 
as opposed to the "negative" gravity, generated by the negative cosmological 
constant. But when cosmological constant becomes considerably large (i.e. 
when his larger than 8irGU or 8nGW or, equivalently, when A is larger than 
U or W) then no soliton energy contribution can prevent this "deviation", 
apart from a shrinking in the soliton magnitude and a consequent decrease 
in the soliton mass and particle number. 

Figure Ogives more interesting results. We start from an u near the crit- 
ical value 1/2 holding true in the absence of gravity. The field configuration 
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is a small one, with small energy and charge. When decreasing u (or A sur ) 
gravity becomes more important and the soliton becomes larger and more 
massive. There is an interesting point, common in similar diagrams of other 
kinds of q-stars. Decreasing frequency both energy and charge increase, but 
only until the frequency reaches a certain point, below which the two quan- 
tities decrease. So for a certain charge there are usually two values of the 
energy, the lower one corresponding to high frequency (lower branch of the 
E = E(Q) diagram) and the higher one corresponding to low frequency (up- 
per branch of the E = E(Q) diagram). So, a soliton star with a certain 
charge and a high energy amount can emit the energy excess, falling down to 
the lower branch, changing of course its frequency. We will study this case in 
more details when treating fermion-scalar q-stars, for which there is an ana- 
lytical solution for a special A ^ case and the obtained relations are simpler 
than in the case of q-stars with two scalar fields. This change in the energy 
and charge variation with respect to the frequency happens only when the 
cosmological constant differs from zero. The above analytical solution given 
for the zero cosmological constant case predicts that both energy and charge 
are monotone decreasing functions of the frequency. When v4 sur — > one 
expects that a black hole is formed. But the calculations depicted in figure Q 
interrupted for an A SUT > value, because below that certain value the star 
decays into free particles as the energetically favorable choice. This means 
that before the formation of an horizon, the q-star decays. So, when A ^ 
gravitational collapse is impossible, at least in the usual way of decreasing 
frequency holding true for q-stars in 3 + 1 dimensions, ^H|-[23j- There may 
be a possibility of forming q-type black holes, if one takes into account the 
spatial variation of the scalar field, but this demands a considerably different 
way of solving the equations of motion to the one known so far. The same 
discussion holds true for any sort of q-star. 



A first simple generalization to the above described soliton is a field configu- 
ration with two scalar fields, one N-carrying <f), and, consequently, complex, 
and the other, a, used to constrain the N-carrying field within a certain re- 
gion, generating an appropriate potential, taken to be real for simplicity. The 
Lagrangian in the case under discussion is: 



3 A q-star with two scalar fields 
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and the metric is supposed to have the same form as in the q-star with one 
scalar field, i.e. to be static and spherically symmetric. The potential has 
the general form: 



We paramertize d as: 



U = a( f) 2 a 2 + b(f) 4 + c(o 2 - df . 



d = ^- 3 /2 fl -l/2 6 l/2 



(3.2) 



(3.3) 



where the parameter p has mass dimensions, is nothing more than a way of 
re-expressing d, is of the same order of magnitude as the potential parameters 
and will simplify the rescalings of the energy quantities. We now rescale the 
matter fields 



= (26) 1 V 4/3 0, a = o 1 / 2 (26)- 1/ V /4 o- • 



We also define: 



A = a 1 bc 

and the potential takes the simple form: 



u = u/p 3 , u = 4> 2 a 2 + ^ 4 + ^{a 2 - i; 



1 74 , 

-6 H — i 
2 r 2 



(3.4) 
(3.5) 

(3.6) 



The rescalings in the other quantities (spacetime and frequency) are the same 
as in eq. 12.161 

We use a simple ansatz for the matter fields, namely we regard that within 
the soliton holds: 

4&t) = <p(p)e-** t , a = 0. (3.7) 
Redefining: f = p87rGp, we have: 



V = A 
with 



( d(p x 



V dp 



da\ 



dp J 



d(p^ 2 



dip^ 



fp) «{^)^«W,U, (3.8) 



W = BCj 2 ^ 2 . 
The Euler-Lagrange equation for the complex field is: 



A 



d 2 <p ( 1 ld{A-B)\ d^ 
df 2 \f 2 df J df 



8nGp 



1 6{U - W) 
2 



(3.9) 



(3.10) 
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Regarding the metric and matter fields as smoothly varying within the star, 
we find from the Lagrange equation for the if field: 



<p 2 = Cu 2 B, (3.11) 



U = -(Cj 4 B 2 + A) 

r ; (3.12) 

W = u 4 B 2 . 

Within the surface the metric and matter fields vary rapidly. Repeating 
the discussion (eqs I2.24H2.28]) of the previous section we find the eigenvalue 
equation for the frequency: 

A ^ /4 



\ sur / 

We will choose A = 1/9 for simplicity. Dropping from now on the tildes, the 
Einstein equations take the form: 

_i4^ = W-±A-A. (3.15) 
2rBdr2 2 y J 

The soliton mass is given by the G° component of the Einstein tensor and 
the particle number is: 

Q = 4vr J drruv 2 ^ = 4vr J drru 3 B 3/2 A~ 1/2 . (3.16) 

Outside the soliton we take ip to be everywhere zero and a equal to unity, 
so that both the energy density and Noether current are zero. These values 
are also solutions to the Lagrange equations of the matter fields. 

We will now find some analytical solutions. The case of zero cosmological 
constant is one of them. In this case holds: 

B(r) = 1/A SUT , A{r) = 1 - 2Ar 2 . (3.17) 

The soliton radius R can be found equal to: 

H -{^f L ) =(^H • (318) 

15 
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Figure 6: The total mass of the field configuration as a function of the particle 
number for several values of the cosmological constant for a q-star with two 
scalar fields. Dashed lines depict analytical solutions. 
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Figure 7: The absolute value of the N— carrying field at the center of the 
field configuration as a function of the cosmological constant for a q-star with 
two scalar fields for four values of A SUT , or equivalently, the frequency. For 
flat spacetime the field has a constant value equal to A 1//4 (here 1/9 1//4 ) in 
the soliton interior. 
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Figure 8: The radius of a q-star with two scalar fields as a function of the 
cosmo logical constant for four values of A SUT . 
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Figure 9: The total mass as a function of the cosmological constant for a 
q-star with two scalar fields for four values of A snr . 
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Figure 10: The particle number as a function of the cosmological constant 
for a q-star with two scalar fields for four values of the A snT . 

and the energy and charge 

E = (1 - A sm )n = (1 - ^ 2 /A 1/2 )vr , (3.19) 

Q = 2v / 3tt(1 - Al£) = 2v / 3tt(1 - cu/A 1/4 ) . (3.20) 

We can see that the equation E = Q has no solution, so, the soliton energy- 
can not be equal to the energy of free particles and consequently the star can 
not decay into free particles. We can also find the relation for the energy E r 
stored within a surface of radius r: 

E r = 2\nr 2 , (3.21) 

which verifies relation l2.25[ when combined with eq. 13.171 

We will now examine a different case admitting analytical solutions. The 

case: 

^A + A = 0. (3.22) 

If R is the soliton radius (i.e. the solution to the equation: A{r) = l/B(r)) 
the soliton parameters take the form: 

R = 3V2(A 2 S £ - A sur ) 1/2 = 3V2(-3u, 2 + 3 2 / 3 ^ 3 ) 1/2 , (3.23) 
(18A SUV -r 2 + R 2 ) 3 (18A1 3 - r 2 ) 3 (-r 2 + 18 • 3 2 / 3 a// 3 ) 3 



A(r) 



5832A2 5832,42 52488w 4 



sur 



(3.24) 
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E = (1 + - 2A sur )rr = (1 - 6u; 2 + 3 2 /^ 4 / 3 )tt 



(3.26) 



Q = 4^(4^ - Al£) = 4^/3tt(3 1 /V/3 - 3^) . 



(3.27) 



The equation E = Q admits analytical solution numerically equal to A sm = 
1.03484 ■ 10 -5 . Then, E = Q = tt and below this value of the surface metric 
the energy of the free particles with the same charge is less than the soliton 
energy, making in this way stars with A sm below that critical value unstable. 
We can also find E r : 



It is a matter of simple algebra to verify equation 12.251 

As one can see from figure HI the analytical solutions behave in the same 
way as the solutions obtained in numerical methods. Summarizing our re- 
sults, we see that decay into free particles is forbidden for the flat spacetime 
case, energetically favorable for an anti de Sitter spacetime, but only below 
a certain value of the frequency. Solitons can not be extremely large. For 
any value of the independent parameters (uj and A) there is a certain region 
in the E-Q phase space, fully depicted in our figures, IBIITUl 

4 Non-abelian q-stars 

The Lagrangian in this case is: 



E r 



2 (792Al 3 + 324A 2 ur - 54A1V + r 4 ) 




(3.28) 




(4.1) 



The Enstein equations are: 



2p dp 



1 dA 




) 



2 



+ A 



(4.2) 



1 AdB 
2pB~dp 




(4.3) 
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with a general renormalizable potential 

U = ^ + |0 3 + ^ , (4.4) 
with in the 5*0(3) 5 representation. We make the following rescalings: 

9 = , = W#)0 , A = # 2 A , p = p> _1 , u = fiw. (4.5) 

All energy quantities are rescaled by p 4 /^ 2 . The potential takes the simple 
form: 

V= p, £f = f + | + ^, («) 
where \ = \/g 2 Redefining 



8nG^rA = A , 

g 2 



the independent Einstein equations take the simple form: 

1 dA 



^ W + W + V + A, (4.7) 

_L^W + V- M -A, (4.8) 



where: 



W = TrC/ , 

W = Tr^£?(^P| =Trfi[fi,0] 2 , 

' 1 «9t y (4.9) 

It is obvious that for a large soliton with smoothly varying fields within the 
interior the V terms are of 0(G) order and, thus, can be neglected. We 
define Vt through: 

=«[«.*]• ( 4 - 10 ) 
20 



We also define: 



1 

n = fi o = -vjfi ( o o o 



(4.11] 



-10 

With the above definitions, the Euler-Lagrange equation takes the form: 



dU 



dcj)J 



d 2 (j) d(j> (I l_dA _ 1 dA\ 
df 2 dr \r 2A df 2A df J 



(4.12) 



e li Vdiage-^. 



where 1 is the unity matrix. We can always diagonalize 
The rigid rotation condition means that R(p, t) = fit + C where we can 
eliminate the constant C through a global 50(3) rotation and write: 

= --^2-diag(l + y,-2,l-y)«»^ = --0 2 .diag(l+j/,-2,l-y) . (4.13) 
The charge is defined by the equation: 

Q = f <PpV=gf = -i j d 2 p^B[4>A] . (4.14) 

The equation of motion for the matter field within the soliton dropping 
the 0(MGp 2 /g 2 ) takes the form: 



+-\<f> — lTr ( -<p z + - 



6 



l ~2 , Itj, 



6 



2Co 2 B( 



} )diag(l,0,-l). (4.15) 



At the inner edge of the surface we obtain the usual relation (repeating eqs 

(^-W) sur = 0. (4.16) 

The usual relation 5(W = U)/5<fi = is no longer valid. We can see that that 
neglecting the 0(8irGfi 2 /g 2 ) in the Euler-Lagrange equation, multiplying by 
4> and tracing yields 2W = Tr(0 ■ dU / d(p). Applying equation 14. 161 we obtain 
for the fields at the inner edge of the surface, with R the soliton radius: 

fa(R) = ™( y2 ~\ ) , (4.17) 

A \(y 2 + 3) / sur 

9 / s fl / 3 \ 3 1 (w 2 - 1) 2 1 
~ 2D ^ m — ' . (4.18) 



3 l (y 2 - iy 

4A y 2 (y 2 + 3f 
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Taking the (2,2) component of the Euler-Lagrange, eq. 14. 151 we take: 



y 2 (f) 



3 \(p 2 2 + 202 + 8 

02 2 — A0 2 



(4.19) 



Using the above equation combined with the 2W = Tr(0 ■ dU/dcj)) relation, 
we find for the 02 field: 



02(0 



A\u 2 B{?) 



2A V L 



4:\ti 2 B(f)\ + 8X[4u 2 B(r) - 1] 



1/2 



Using now eqs 14.171 14.191 we find: 

[(A - \){y 2 + 3) 3 + 16(y 2 + 3) 2 - 72(y 2 + 3) + 







(4.20) 



(4.21) 



A < 1 holds in order 14.191 to have real solutions. The solution to the above 
equation is: 

v\R) = 16 



2 2 1 / 3 
3-1 + A 



2 2 1 / 3 

+ 77- 



r,2 



-13 + 162A-81A 2 + 9(-l + A) 



\ 



3(1 -A) 

1/3 



-1 + 81A 2 
(-1 + A) 2 



31- A 



Ts2 



13 + 162A + 81A 2 + 9(-l + A) 



81A 2 



\ (-1 + A) 2 



+ 

1/3 



(4.22) 



The idea for the non-abelian solitons is the following: We give a certain value 
to the free parameter A, in our calculations A = 2/3, and using eq. 14.221 we 

find y 2 (R)- Substituting in eq. 14. 181 we find the eigenvalue for the frequency, 

1/2 ~ 
in our case equal to 0.4956A S ur- We then substitute the values of the 02 

and y2 matter fields (eqs l4.19H4.20j) in Einstein equations and solve them 

numerically. 

From the figures fTTlfTol we see that for larger values of the cosmological 
constant the matter field (02) is larger, but the soliton radius, charge and 
energy decrease with the increase of the cosmological constant in absolute 
values as in other kinds of q-stars. 

We now investigate the analytical solution to the Einstein equations in 
the case of zero cosmological constant. We find: 



.4 



2.82529f 2 



(4.23) 
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Figure 11: Energy as a function of the particle number for four values of the 
cosmo logical constant, for a non-abelian q-star. 
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Figure 12: The matter field value (02(0)) at the center of the field configu- 
ration as a function of the cosmological constant for a non-abelian q-star for 
four values of the A SUT . 
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Figure 13: The soliton radius as a function of the cosmological constant for 
a non-abelian q-star for four values of the A SUT . 
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Figure 14: The total mass as a function of the cosmological constant for a 
non-abelian q-star for four values of the A SUT . 
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Figure 15: The particle number as a function of the cosmological constant 
for a non-abelian q-star for four values of the A sm . 



R=[(l- A sur )/2.8259] 1/2 , (4.24) 
E = n(l-A SUI ) , (4.25) 
Q = 16.7482(1 - Alii) ■ ( 4 - 2 6) 

5 Fermion- Scalar q-stars 

Fermion-scalar q-stars are realistic field configurations of a scalar and a 
fermionic field. These field configurations are supposed, [21], to describe 
certain stellar objects known as neutron stars. The fermionic field carries 
the charge needed to stabilize the soliton and the real scalar field generates 
an everywhere positive potential. In the star interior the scalar field is ap- 
proximately zero but outside the soliton takes a certain value o"o, vanishing 
the potential. The potential energy density outside the soliton vanishes in 
this way (eq. 15. 9j) and the fermion mass acquires its vacuum value m. 

We regard the fermion-scalar q-star as a zero temperature fermionic sea 
with local Fermi energy and momentum e 2 F + k F + m 2 (cr). With spherical 
symmetry, the local fermion density is: 

m = i L/ n h *k = g , (5.1) 
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and the scalar density: 

, - , 1 f , Q , m m(—m + EF) 

M) = — / n k d 2 k = — . 5.2 

We also find using the definitions for the energy and pressure density: 

If k 2 1 

P^ = - n k d 2 k— — ^ = ~(e F {^ip) - miipip)) , (5.3) 

z J 2(Ar + m z ) ' <J 

^ = -L y n fe rf 2 A;(A; 2 + m 2 ) 1/2 = 2P^ + m(#) . (5.4) 
The Lagrangian density is: 

£ = - m{a))^ - i(^) 2 - U(a) , (5.5) 

with 

m(a) = ga . (5.6) 

We choose g = 2 for our calculations. The Euler-Lagrange equation for the 
scalar field is: 

d 2 a Ida <9 . __ N _ N 

^ + ^ = -a; (p *- c/) - (57) 

An appropriate form for the potential is: 

17 = \\{a 2 - a 2 f. (5.8) 

The proper generalization in order to include gravity in our discussion is 
to replace the Fermi energy with a global chemical potential, according to 
the equation I5.14[ [T2*] . In this way the Fermi energy takes the proper gen- 
erally covariant form, being the zero component of the energy-momentum 
4— vector. The generalization of the Lagrangian with the inclusion of gravity 
is straightforward: 

£/V=9 = \{4>1^ - VVT'V) - m{<r)W + \a-,^ - U(a) . (5.9) 

The proper rescalings for the case under discussion are made with respect 
to a new mass scale \x equal to a 2 , simplifying in this way the potential 
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and extracting a fi 3 factor from the Lagrangian. So the quantities in the 
Lagrangian are rescaled in the following way: 

a = — , ^ = — , A = — , m = — , p = pa , £ = — . (5.10) 

CT cr cr tr CT 

We also make two new redefinitions: 

r = 8nGa 2 , A = -A, . (5.11) 

07TG(Jq 

So, dropping the tildes we can regard that inside the soliton the matter 

field derivative with respect to the radius is of order of 0(Gfi 2 ) so the Euler- 
Lagrange equation for the scalar field can be written: 

^-(P^-U) = 0. (5.12) 

The energy-momentum tensor is: 

= -t($yby) - ^yty) + ° 4l o ,v ~ 9" u [l/2a. a a' a - U{a)\ . (5.13) 

We define a global chemical potential through the equation: 



l = elB-\r) . (5.14) 



This quantity retains a role similar to the frequency of the other kinds of 
soliton stars. The particle number is: 

Qip = J ^^dx 1 dx 2 j° = 2tt J rdry^(^V) • (5-15) 

The Einstein equations take the form: 

1 dA , . 

-— — = £^ + U + A, 5.16 

2r ar 

~hw§ : = p *~ u ~ k ' ^ 

Within the surface we repeat the discussion of eqs. 12 .24112. 251 and find 
from eq. 15.121 
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Figure 16: Energy as a function of charge for four values of the A SUT or, 
equivalently, the chemical potential for a fermionic q-star. We start from a 
small value of the cosmo logical constant. Both energy and charge have small 
values. When increasing the cosmological constant, in absolute values, both 
energy and charge increase in order the soliton to generate "positive" gravity 
as a counter effect to the "negative" gravity from the spacetime structure. 
But when A takes very large absolute values, a large soliton can not be stable 
and mass, charge and radius decrease as one can see from figures ITMT^l 



This is the eigenvalue equation for the chemical potential. The Einstein 
equations take the soluble form: 

1 dA 1 q _ Q/n A 



2r dr 6i * I 
1 AdB 1 



^,£? 3 / 2 + - + A , (5.19) 
r.l^ - 7 - A , (5.20) 



2r B dr 12vr 
The total charge is: 

Q = \J drrulBA- 1 ' 2 . (5.21) 

For our numerical calculations we use A = 1/2. 

Outside the soliton we regard that a is everywhere unity and ip zero, so 
both energy density and current are zero. 

Let us now take A — > 0. In that case the Einstein equations admit the 
following analytical solution, also verified numerically: 

/ Q \ 2/3 q 

B{r) = l/A snr ={-\ u-\ A(r) = l--r\ (5.22) 
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Figure 17: The soliton radius as a function of the chemical potential for a 
fermionic q-star for five values of the cosmological constant. When toy, is 
large the gravity strength according to the eq. 15. 181 is small. 
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Figure 18: The total soliton mass as a function of the chemical potential for 
five values of the cosmological constant for a fermionic q-star. 
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Figure 19: The particle number as a function of the chemical potential for 
five values of the cosmological constant for a fermionic q-star. 




Figure 20: The total mass of the field configuration as a function of the charge 
for a fermionic q-star with constant cosmological constant when varying the 
chemical potential. We start from large values of this quantity (i.e: with A SUT 
close to unity, according to eq. 15. 18 J) and, consequently, weak gravity Then 
both charge and energy have small values. Decreasing both energy and 
charge increase up to a certain value of the chemical potential, below which 
they both decrease. 
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We can easily find that the parameters of the soliton are: 



R=[(l- A sur )8/3] 



1/2 



2 \ 



2/3 



1/2 



£ = (1 - A sur )7T = 



2 \ 2/3 

3tt/ 



7T 



3vr\ 2/3 4 



3vr\ 2/3 4 



2 ) 3 



/2\ 1/3 



(5.23) 



(5.24) 



(5.25) 



where R is the soliton radius, E the total mass and Q the soliton charge. 

We will now examine a different case admitting analytical solutions, 
namely: 

A = — - (5.26) 
4 

In that case Einstein equations take the form 



1 dA _ 2 3p3/2 

2r rfr 12tt ^ 



1 



2r 5 rfr 12tt" ^ 
We can find the soliton radius: 



£3/2 



R = ^-A SUI ) 1/2 

The metrics take the form: 



2\* 



A(r) 



(16A sm - r 2 + i? 2 ) 



65536AL 



(16Al 4 - r 2 ) 4 _ vr 2 (3r 2 - 16^^) 



65536AL 



2359296^ 



(5.27) 
(5.28) 

(5.29) 



(5.30) 



B(r) 



1 (16A sur - r 2 + R 2 f 
256 " " .41,. 



1 (16A 3 u 7 r 4 -r 2 ) 2 (37rr 2 - WV^uff 



256 



A3 
sur 



1024c^ 



(5.31) 
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We can also find the energy and charge: 

E = (1 + 2A^ - 3A sur )n = n - 2 2 / 3 (3vr) 1 / 3 ^ + 2^| Jf , (5.32) 

Q = 2^(3nf \AHt - Ag) = 2^uf - 2 5 / 3 (3 7 r) 1 % . (5.33) 

For both the analytically soluble cases the chemical potential varies from 
the maximum value obtained by the eq. 15.181 (equal to 37r/2 for these stars) 
when the minimum value is obtained equating the soliton energy with the 
energy of the free particles with the same charge. When B = A~^ = const, 
this equation has no solution and no decay to free particles is possible. In 
the latter case the minimum value is u)^ = 0.0481. 

We will now examine the A + A/4 = analytical solution in more details 
because in the fermion-scalar soliton stars we face the simpler relations when 
compared with those resulting from the investigation of the solitons with 
two scalar fields. With the help of figure EH we see that the E = E(Q) 
function has two branches, a lower corresponding to large values of A sm , 
or, equivalently, u^, and the upper one corresponding to small values of the 
above quantities. The boundary separating the two areas is determined by 
the solution of the equations 

The solution to both of the above equations is only one: 

< = j(f) 1/3 , ^ r = i. (5-34) 

For uj^ > cuT energy and charge are monotone decreasing functions of the 
chemical potential, but they are both monotone increasing for uo^ < u™. The 
charge shows a more rapid decrease. 

Solving eq. I5.33l in terms of A SUT for simplicity we find two solutions, called 
Asii and A^ corresponding to the upper and lower branches respectively of 
the the E = E(Q) function. These solutions are: 

A W = J_ ( 21/332/3^5/3 _ 12 . 2 2 / 3 3 1 / 3 Qtt 7 / 3 + 36tt 3 - 6V2tt 2 x 

\/-2Q 3 - 12 • 2 2 / 3 3 1 /3Q vr 4/3 + 18tt 2 + 5 • 2 1 / 3 Q 2 (3n) 2/3 J (5.35) 
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A® = J— ( 2 1 / 3 3 2 / 3 Q 2 vr 5 / 3 + 12 • 2 2 / 3 3 1 / 3 Qvr 7 / 3 + 36tt 3 - 6v / 2vr 2 x 
727r 3 V 



-2Q 3 - 12 ■ 2 2 / 3 3 1 /3Q vr 4/3 + lg7r 2 + 5 . 2V3g2( 37r )2/3j ^ 36 ^ 

The metric of the eq. 15.351 is the small one (less than 1/16) when the the 
metric of eq. 15.361 is the large one. Substituting in eq. 15.321 we find the 
energies corresponding to the two branches: 



2 2 ; 4-23/2(3vr) 1/3 



-yJ-AQ 3 - 24 • 2 2 / 3 3 1 /3g vr 4/3 + 367r62 + 1 . 2 1 ' 3 Q 2 (3 



nf 3 + 



' (2 1 / 3 3 2 / 3 g 2 vr 5 / 3 - 12 • 2 2 / 3 3 1 / 3 gvr 7 / 3 + 6tt 2 x 



6 ■ 2 1 /4 3 l/2 7r 5/4 

\ 3/4 

vr _ J-AQ 3 - 24 ■ 2 2 / 3 3 1 /3Q vr 4/3 + 36tt 2 + 10 • 2 1 /3Q2( 37r ) 2 /3) ^5 37) 



the large value of the energy corresponding to the upper branch and: 

2 +C ^2; 4-2 3 / 2 (3vr) 1 / 3 
Iy / _4g 3 _ 24 • 2 2 / 3 3 1 /3Q 7r 4/3 + 367r62 + 10 . 2 1 / 3 Q 2 (3n) 2/3 + 
1 



7T 2 X 



V'^QV" - 12 • 2^3'/ V/3 + 6, 

I \ 3/4 

(6tt + \J-AQ 3 - 24 • 2 2 / 3 3 1 / 3 Q 7 r 4 /3 + 36vr 2 + 10 • 2 1 / 3 Q 2 (3tt) 2/3 ) J ,(5.38) 
the small one. 

6 Conclusions 

Our discussion can be summarized as follows: 

1. There are stable, gravitating field configurations, in 2 + 1 dimensions, 
consisted of a iV— carrying field, fermionic, or scalar (abelian or non-abelian) . 
An always positive potential, generated either by the the N— carrying field, 
or by an additional real scalar field, is necessary in order to stabilize the 
soliton. The potential as a function of the field increases as <j) 2 for small 
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values of the field, then slower than 2 providing in this way a local minimum 
for some value of the field differing from zero and for large values of the field 
as 4 , or 6 usually, in order to ensure the positivity of the potential for large 
values of the field. The soliton solution lies near this minimum and the total 
energy of the star is smaller than m 2 2 , i.e. the potential energy of the free 
particles in any field configuration under consideration . There is always a 
certain frequency, minimizing the total energy of the configuration, obtained 
by the equality of the W— energy to the [/—energy This equality holds true 
in the absence of gravity and is generalized in a generally covariant way. 

2. When A 7^ and A SUT — > the particle number decreases rapidly but 
the total energy decreases slowly, as one can see from figures [H |Hl EU and 
I2TH So, the A snr = case (formation of horizon) is impossible, because below 
a certain value of A sm fission into free particles is energetically favorable. 
This result is supported by both numerical and analytical solutions. Decay 
into free particles is energetically forbidden for zero cosmological constant, 
as can be proved by the analytical solutions for every special kind of q-star. 
For cosmological constant differing from zero, the energy and the charge of 
the soliton increase when the gravity becomes stronger (i.e.: when A SUT de- 
creases) but below a certain value of the frequency energy decreases slowly 
and charge more rapidly. The result is that for a certain charge there are 
usually two different values of the energy. So the soliton with the larger en- 
ergy can emit the energy excess and goes to the lower-energy state. Due to 
the above behavior of the E = E(Q) function, we found analytically and nu- 
merically that below a certain value of the frequency the soliton decay to free 
particles is energetically favorable, so when the cosmological constant differs 
from zero there is a certain range of the soliton parameters (the independent 
parameter, frequency, and the derivative ones, radius, energy and charge) 
which verify the soliton stability. This range of the phase space parameters 
is fully depicted in our figures. 

3. When the cosmological constant is zero or small with respect to the 
energy densities (namely U and W), the soliton parameters, radius, mass 
and particle number, increase when the frequency decreases, i.e. when the 
gravity becomes stronger. When the cosmological constant becomes more 
important than the energy densities, the soliton becomes smaller and the 
mass and particle number show a corresponding decrease. 
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